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Abstract 

We study the scintillation index of partially overlapping collimated lowest order Gaussian laser 
beams with different wavelengths in weak atmospheric turbulence. Using the Rytov approximation 
we calculate the initial beam separation that minimizes the longitudinal scintillation. Further 
reduction of the longitudinal scintillation is obtained by optimizing with respect to both beam 
separation and spot size. The longitudinal scintillation of the optimal A^-beam configurations is 
inversely proportional to N, resulting in a 92% reduction for a 9-beam system compared with 
the single beam value. The radial scintillation values for the optimal A^-beam configurations are 
significantly smaller than the corresponding single beam values. 
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I. INTRODUCTION 



Propagation of light through atmospheric turbulence is the subject of intensive research 
owing to the many applications in free space laser communications l|. In these applications 
it is desirable to reduce the turbulence effects on the propagating light. One promising pos- 
sibility to achieve this goal is by using temporally partially coherent optical field consisting 
of multiple laser beams with different wavelengths 0, ^, Indeed, in a typical setup in 

which the detector's response time is large compared with the inverse of the frequency dif- 
ference between any pair of beams in the input field, rapidly oscillating contributions to the 
total intensity average out. Consequently, one can expect smaller values of high moments of 
the intensity compared with corresponding single beam values. This would result in smaller 
values for the scintillation index and for the average signal to noise ratio (SNR). 

Generation of temporally partially coherent light consisting of multiple beams with dif- 
ferent wavelengths can be efficiently realized by using an array of vertical external cavity 
surface lasers (VECSELs). These devices have the advantage of generating high power, 
spectrally narrow, wavelength tunable TEMqo beams (lowest order Gaussian beams) ^. 

Propagation of temporally partially coherent light in atmospheric turbulence was first 
studied by Fante [tI, Q] , who obtained approximate analytic expressions for the scintillation 
index of a single infinite planar wave. More recently Kiasaleh studied propagation of an 
infinite multi-wavelength planar wave in weak atmospheric turbulence and showed that the 
achievable SNR is larger in the multi-wavelength case compared with the single-wavelength 
case [21, |3|]. These previous studies focused on infinite planar waves, whereas in practice, 
Gaussian laser beams with finite initial spot size are employed. Since the dynamics of 
the optical field can strongly depend on the initial spot size it is important to take into 
account the finite spatial dimension of the beams. Furthermore, the optical field of 
collimated TEMqo beams with different wavelengths depends on the wavelength separation 
and also on the spatial separation between the beam centers at the transmitter. Therefore, 
one can exploit these two different dependences to reduce scintillation and optimize system 
performance against turbulence effects. The dependence of the scintillation index on the 
wavelength separation of multiple overlapping TEMqo beams was studied in Ref. [4]. It was 
found that for typical lasercom setups modest scintillation reduction of about 10% can be 
acheived by controlling the wavelength difference between the beams. We emphasize that 
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only the case where the beams are completely overlapping at the transmitter plane was 
considered in Ref. ^|. In the current Letter we focus attention on scintillation reduction 
by varying the spatial separation between the beams at the transmitter. We show that this 
approach leads to a much stronger decrease of the scintillation compared with the approach 
employed in Ref. ^. Moreover, our approach allows us to find the initial beam separation 
that minimizes scintillation, thus providing a simple solution for the important problem of 
optimizing temporally partially coherent sources of light against turbulence effects. 

In Ref. jsl we took the first step in this approach and established the framework for 
calculating the scintillation index for multiple partially overlapping beams in weak atmo- 
spheric turbulence. Using the Rytov approximation and considering a typical 2-beam sys- 
tem we found the initial beam separation that minimizes the longitudinal scintillation. We 
showed that the longitudinal and radial scintillation for the optimal 2-beam configuration 
are smaller by about 50% and 35%-40%, respectively, compared with the corresponding 
single-beam values. However, two important aspects of the problem were not addressed in 
Ref. jsl]: (1) the A^-dependence of the longitudinal scintillation reduction compared with the 
single-beam case, and (2) the possibility to optimize the system with respect to both initial 
beam separation and initial spot size. In this Letter we address these two central issues in 
detail. 



II. CALCULATION OF THE SCINTILLATION INDEX 

Consider propagation of collimated linearly polarized TEMqo beams with differ- 
ent wavelengths Xj, j = 1,...,N, in weak atmospheric turbulence. Assuming that the 
beams propagate along the z axis and denoting by the beam-center locations at the 
input plane z = 0, the magnitude of the total electric field E at z = is E{r,0,t) = 
'^f=i Uj{Tj, 0) exp [—iujjt], where Uj{rj, 0) = exp[— r|/14^Qj], r is the radius vector in the xy 
plane, rj = r — dj, t is time, kj = 27[/Xj are wavenumbers, ujj = kjc are angular frequencies, 
and c is the speed of light. In addition, W^j are the initial spot sizes and we assume that all 
beams have the same amplitude. Assuming weak turbulence, the propagation is described 
by N uncoupled linear wave equations 

\7^Uj + [1 + 2ni(r, z)] = 0, (1) 
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where ni(r, represents the refractive index fluctuations, |ni(r,2;)| ^ 1. To solve Eq. ([T]) 
we follow Ref. [l| and employ the paraxial approximation together with the Huygens-Fresnel 
integral and the second order (with respect to rii) Rytov perturbation method. The total 
intensity at 2; = L is 

N 



N N 

5^^f/,(r,,L)[/;(r^,L)exp 



(2) 



where Ij{rj,L) = \Uj{rj, L)\'^ is the intensity of the j-th beam. The intensity measured by 
the detector is the time average /dei(r, L) = r'^ dt/(r, L, t), where r is the response time 
of the detector. Assuming a slow detector and 7^ for j 7^ m, we neglect the terms 
UjU^ j 7^ m, which are rapidly oscillating with time. Therefore, the measured intensity is 

N 

/,e*(r,L)^ J]/,(r„L). (3) 



The total scintillation index for the A^-beam system is 

a?(r,L) = (/L(r,L))/(/,et(r,L))2 



(4) 



where (...) stands for average over different realizations of turbulence disorder. Using Eqs. 
([3]) and (jlj) we obtain 
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(5) 



The total scintillation index a] can be decomposed into a longitudinal component aji{L) = 
erf (0, L) and a radial component cr^(r, L) = af (r, L) — aj i{L). 

In calculating intensity moments we assume that the perturbation field in the Rytov 
approximation is a Gaussian random variable and that the turbulence is statistically ho- 
mogeneous and isotropic. Consequently, the average intensity of the ?'-th beam is given by 





2r| 



Hij{rj,L) 



(6) 



where Wj is the spot size at distance L, and Hij is expressed in terms of a double integral 
of the spectral density of the refractive index fluctuations over wavenumber k and 

propagation distance z. [See Ref. js], Eq. (18)]. The average of the second moment is {ll 

(/J(r„ L)) = (/, (r„ L)f exp [H^.ir,, L)] , (7) 

where is given by another double integral of ^n{i^) over k and z. [See Ref. jSj, Eq. 
(20)]. The cross-intensity term {Ij{rj, L)Im{'i^m, L)) is given by 

{Ij{rj,L)Im{rm,L)) = (/^(rj, L))(/m(r„, L)) x 

exp {i?2jm,(rj, Tfn] kj, km) + E2mj(jm, ^j] ^mj ^i)"l" 

2Re[E3jm{rj,rm;kj,km)]} , (8) 



where -E'2jm, -E'2mj and E^jm are three different integrals of $„(k) over k and z. [See Ref. 
Is], Eqs. (22-24)]. 

We consider two typical free space laser communication setups, in which the central 
wavelength is Ac = lO^^m, the wavelength spacing is AA = lO^^m, all beams are collimated 
and have the same initial spot size and on-axis amplitude. In setup A L = 1km and the 
refractive index structure parameter is = 3.0 x 10~^^m~^/^, and in setup B L = 10km 
and = lO^^^m"^/'^. Both setups correspond to weak atmospheric turbulence conditions, 
where the Rytov variance aj^ = 1.23C^fc^/6L"/6 is about 0.1 and 0.23, respectively, for all 
beams. We use the Von Karman spectrum to describe the refractive index fluctuations. 
Thus, $„(/€) = 0.033C2 (/t2 + /t2^J~"/%xp {-k^/k^J, where = 5.92//o, Hout = l/^o, 
and Lq are the turbulence inner and outer scales, respectively, and Iq = 1.0mm, Lo=1.0m 
are used. 

For even we consider initial configurations in which the beam centers are located on 
a circle with diameter d centered about the ^-axis with equal angles between dj_i and dj. 
Thus, for = 4, for example, the centers are at di = dSt/2, di = d'y/2, di = —dSt/2, and 
di = —d'y/2. For odd we consider the same geometry as in the N — 1 case, with an 
additional beam on the z-axis. 

The (i-dependence of the longitudinal scintillation for 2-, 3-, 4-, 5-, 8-, and 9-beam con- 
figurations in setup A with initial spot sizes Wq = 1cm is shown in Fig. [1] together with 
the corresponding value for a single TEMqo beam with the same total power and initial 
spot size. One can see that in each of the N > 2 cases the crf^-curve exhibits a minimum 
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FIG. 1: Longitudinal scintillation cr| ^ vs initial beam separation d for setup A. The solid, dashed, 
dotted, dashed-dotted, short-dashed, and short-dotted lines correspond to 2-, 3-, 4-, 5-, 8-, and 
9-beam configurations, respectively. The square stands for the value for a single beam with the 
same total power and initial spot size. 

at an intermediate d value, do = 2.8cm, 4.4cm 3.6cm, 4.6cm, 4.0cm and 4.8cm, for the 
2-, 3-, 4-, 5-, 8-, and 9-beam configurations, respectively. These minima correspond to the 
optimal configurations of the N beams for the given physical conditions and geometric ar- 
rangements, where optimization is with respect to longitudinal scintillation. Comparison 
with the single-beam value shows that the longitudinal scintillation is reduced by 53.4%, 
63.2%, 82.2%, 84.4%, 88.1%, and 92.0% for the 2-, 3-, 4-, 5-, 8-, and 9-beam optimal con- 
figurations, respectively. Moreover, analysis of the longitudinal scintillation values for the 
optimal configurations shows that a] i decreases like with increasing A^. 

An important question concerns the possibility to further reduce the scintillation by 
optimizing with respect to the initial spot size. This question is addressed in Fig. [21 which 
shows the d-dependence of the longitudinal scintillation for different Wq values for a 5-beam 
system in setup B. One can see that the minimum value of a] i first decreases with increasing 
Wq and then increases. Hence, the optimal configuration for the 5-beam system in setup B, 
is the one with Wq ~4cm and do = 13.2cm. Notice that the final free space spot size of the 
beams for the optimal configuration in this case is only 8.8cm. 

When the spot size is comparable with the radius of the receiver's collecting lens, the 
radial scintillation becomes important. In this case it is essential to understand whether the 
radial scintillation of the optimal A^-beam configurations, where optimization is with respect 
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FIG. 2: Longitudinal scintillation ^ vs initial beam separation d for different Wq values for a 5- 
beam system in setup B. The solid, dashed, dotted, dashed-dotted, short-dashed, and short-dotted 
lines correspond to Wo=0.5cm, 1cm, 2cm, 4cm, 5cm, and 6cm, respectively. 
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FIG. 3: Circularly averaged radial scintillation index cj^^ vs radius r for the optimal 4-beam 
(squares) and 2-beam (circles) configurations, and for a single beam with the same power and 
initial spot size (triangles), in setup A. 

to longitudinal scintiWation, is sufficiently small compared with the single-beam value. Notice 
that in the A^-beam case af.{r, L) is not radially symmetric. To enable comparison with the 
single-beam case we define the circularly averaged radial scintillation in the A^-beam case as 
cr^^(r, L) = ((^^(r, L))g, where {■ ■ .)e denotes averaging over the angle 9. Figure [3] shows the r- 
dependence of af.^. for the optimal 4- and 2-beam configurations and for a single beam with the 
same power and initial spot size, all in setup A. One can see that the radial scintillation for 
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the optimal 4-beam and 2-beam configurations are smaller by about 65%-80% and 35%-40%, 
respectively, compared with the corresponding single beam values. Therefore, optimization 
of the A'"-beam configurations with respect to the longitudinal scintillation leads to significant 
reduction in the radial scintillation, and this reduction effect grows with increasing A^. 

III. CONCLUSION 

We calculated the scintillation index for an array of N partially overlapping coUimated 
TEMoo beams with different wavelengths in weak atmospheric turbulence using the Rytov 
perturbation method. We showed that both the longitudinal and the radial scintillation can 
be significantly reduced compared with the corresponding single-beam values by optimizing 
the beam array with respect to initial beam separation and spot size. These reduction efi^ects 
grow with increasing N, resulting in a 92% reduction in the longitudinal scintillation for an 
optimal 9-beam system. 
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rial Command, USAF, under grant AFOSR FA9550-04- 1-0213. JVM acknowledges financial 
support from the Alexander von Humboldt Foundation. 



[1] L. C. Andrews and R. L. Phillips, Laser Beam Propagation through Random Media (SPIE 

Press, Bellingham, Washington, 1998). 
[2] K. Kiasalch, J. Opt. Soc. Am. A 21, 1452 (2004). 

[3] K. Kiasaleh, "Impact of turbulence on multi-wavelength coherent optical communications," in 
Free-Space Laser Communications V, D. G. Voelz and J. C. Ricklin, eds., Proc. SPIE 5892, 
58920R1 (2005). 

[4] K. Kiasaleh, J. Opt. Soc. Am. A 23, 557 (2006). 

[5] A. Pclcg and J. V. Moloney, J. Opt. Soc. Am. A 23, 3114 (2006). 

[6] L. Fan, M. Fallahi, J. T. Murray, R. Bedford, Y. Kaneda, A. R. Zakharian, J. Hader, J. V. 

Moloney, W. Stolz, and S. W. Koch, App. Phys. Lett. 88, 0211051 (2006). 
[7] R. L. Fante, Radio Sci. 12, 223 (1977). 
[8] R. L. Fante, J. Opt. Soc. Am. 69, 71 (1979). 



8 



